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ABSTRACT 

7/  -  V  . 


A  new  methodology*1ias  been  developed  to  predict  the 
impulsive  noise  generated  by  a  transonic  rotor  blade.  The 
formulation  uses  a  full-potential  finite-difference  method  to 
obtain  the  pressure  field  close  to  the  blade.  A  Kirchhoff 
integral  formulation  is  then  used  to  extend  these  finite- 
difference  results  into  the  far-field.  This  Kirchhoff  foumuU, 
,  t»  written  in  a  blade-fixed  coordinate  system,  ^requires 
initial  data  across  a  plane  at  the  sonic  radius.  This  data 
is  provided  by  the  finite-difference  solution.  Acoustic  pres¬ 
sure  predictions  show  excellent  agreement  with  hover  exper¬ 
imental  data  for  two  hover  cases  of  0.88  and  0.90  tip  Mach 
number.  The  latter  of  which  has  delocalised  transonic  flow. 
These  results  represent  the  first  successful  prediction  tech¬ 
nique  for  peak  pressure  amplitudes  using  a  computational 
code. 
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LIST  OF  SYMBOLS 

<iq  =  atmospheric  speed  of  sound,  m/s 
e  =  rotor  chord  length,  m 

h  =  nondimensional  time  step 

/  =  identity  matrix 

J  =  Jacobian  of  coordinate  transform  matrix 

,Vf  =  hover  tip  Mach  number 

Mi  =  local  Mach  number 

p  =  absolute  pressure,  N/m* 
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P  =  nondimensional  scaled  disturbance  pressure, 

(j *  -  Po)/(*»o«o«*/S) 

po  =  atmospheric  pressure,  N/m1 

r  —  dimensional  distance  along  the  rotor  in  the 

s panwise  direction,  m 

R  =  radial  distance  to  the  rotor  tip,  m 

Res  =  free  stream  subtraction  term  in  Eqn.  (5) 
t  =  nondimensional  time  in  nontransformed  coordi¬ 

nate  system 

t'  =  dimensional  time,  s 

T  =  scaled  asimuthal  angle  location, 

(1/0(0  +  0  —  tan-1  0) 

U  =  contravariant  velocity  normal  to  the  q,  f  plane 

V  =  contravariant  velocity  normal  to  the  £,  (  plane 

Vr  =  nondimensional  scaled  radial  velocity  compo¬ 

nent,  jr/wRt 

W  =  contravariant  velocity  normal  to  the  £,q  plane 
x,  y,  s  =  nondimensional  coordinate  system  aligned  with 
the  blade 

0  =  blade  radial  location,  ~  Ij1^1 

4  =  f-'/J 

i  =  inverse  of  the  blade  aspect  ratio,  e/ R 

1  =  specific  heat  ratio 

T  =  jump  in  potential  across  the  wake  normalized 

by  («oc) 

u  =  blade  angular  speed,  rad/a 

$  =  dimensional  velocity  potential,  m*/s 

♦  =  non-dimensional  velocity  potential,  d/coo 

p  =  fluid  density  normalised  by  free-stream  values 

p  =  fi/J 

p0  =  atmospheric  density,  kg/m 3 
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a  —  nondimensional  integration  variable,  -z 

t  =  normalized  time  in  the  transformed  coordinate 

system 


_  =  azimuthal  angle  in  cylindrical  coordinates  rela- 

;  r  Sivw’to  die  blade 
l  '=  transformed  coordinate  system 
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Delocalization  occurs  for  a  high-speed  rotor  when  the 
supersonic  flow  region  at  the  blade  tip  connects  to  the 
farfidtd'beyond  the  sonic  cylinder.  The  sonic  cylinder  is  the 
radial  location  beyond  which  the  undisturbed  flow  appears 
to  be  supersonic  to  an  observer  on  the  blade.  Flowfield  non- 
linearities  are  very  significant  at  delocalization1  and  these 
nonlinearities  cause  large  changes  in  the  far-field  acoustic 
levels. 


Linear  methods  with  monopoles  and  dipoles  have  con¬ 
sistently  underpredicted  the  far-field  acoustic  levels  for  high 
tip  Mach  numbers.  The  methods  developed  by  Ffowcs 
Williams  and  Hawkings2  improve  upon  linear  formulations 
by  including  the  nonlinear  quadrupole  term  in  the  acoustic 
pressure.  Several  methods  have  been  developed  for  evalu¬ 
ating  this  quadrupole  term3-5,  but  difficulties  concerning 
both  the  size  of  the  quadrupole  volume  and  the  acquisition 
of  aerodynamic  data  make  this  term  difficult  to  integrate 
accurately.  To  date,  neither  computer  codes  nor  experi¬ 
mental  techniques  have  been  able  to  successfully  predict 
the  far-field  acoustic  pressures  associated  with  high-speed, 
delocalized  rotor  flows. 


Computational  fluid  dynamics  (CFD)  computer  codes 
in  their  non-appraximated  forms  properly  handle  all  flow- 
field  non-linearities.  Several  computer  codes  exist  for  solv¬ 
ing  the  three-dimensional  full-potential  equation  for  heli¬ 
copter  rotor  flowfields.9-10  Typically  these  codes  are  well 
developed  for  surface  pressure  and  load  predictions  close 
to  the  blade.  Far  from  the  blade  however,  the  grid  den¬ 
sities  are  insufficient  to  resolve  the  details  of  the  acoustic 
pressure  field.  Extending  these  codes  to  predict  acoustics 
is  too  costly  and  perhaps  impossible  due  to  the  huge  grids 
required  and  instability  problems  that  result  from  increas¬ 
ing  the  size  of  the  computational  domain.  These  instabil¬ 
ities  occur  as  the  radial  distance  away  from  the  blade  is 
increased,  thereby  increasing  the  freestream  Mach  number 
relative  to  the  blade.11 


The  idea  in  this  paper  is  to  use  a  combined  finite- 
difference  code  and  Kirchhoff  formulation  in  order  to  pre¬ 
dict  the  acoustic  far  field.  The  finite-difference  code  used  is 


FPR,  which  stands  for  Full-Potential  Rotor  code.  It  solves 
the  fully  conservative,  three-dimensional,  unsteady  form  of 
the  full-potential  eqauation  and  is  described  in  Refs.  8  and 
9.  The  FPR  code  is  used  to  predict  the  pressure  field  both 
on  the  blade  and  along  a  plane  located  at  the  sonic  radius. 
A  computed  solution  at  the  sonic  radius  is  well  within  the 
limits  imposed  by  the  maximum  grid  size  and  instability 
problems. 

In  extending  this  result  to  the  far  field,  we  use  a 
Kirchhoff  integral  formulation  that  has  been  developed  by 
Isom12.  This  Kirchhoff  formulation  uses  integrals  of  the 
finite-difference  results  at  the  sonic  cylinder  to  determine 
the  acoustic  pressures  in  the  far  field.  It  treats  the  linear 
sonic  cylinder  as  a  caustic  where  all  the  important  data 
originates.  This  data  is  then  propagated  along  linear  char¬ 
acteristics  to  the  far  field. 

The  Kirchhoff  formula  uses  non-rotating  surface  inte¬ 
grals  of  pressure  and  velocity.  This  non-rotating  control 
surface  should  be  luge  enough  to  enclose  the  rotor  and 
the  region  of  nonlinear  flow.  The  linear  sonic  cylinder  suf¬ 
ficiently  encircles  the  nonlinear  tip  region.  The  flow  field 
on  the  sonic  cylinder  comes  from  a  numerical  solution  us¬ 
ing  a  conservative  full  potential  transonic  code.  This  ap¬ 
proach  has  several  practical  advantages.  A  large  volume 
of  quadrupole  integral  is  replaced  by  a  surface  integral 
over  the  sonic  cylinder  surface.  This  control  surface  is  a 
non-rotating  one,  eliminating  mathematical  singularities  of 
the  Ffowcs  Williams  and  Hawkings  equation.  More  impor¬ 
tantly,  the  close  proximity  of  the  sonic  cylinder  to  the  rotor 
tip  means  that  the  finite-difference  code  only  needs  to  solve 
the  near-field  region  near  the  tip. 

THE  FULL-POTENTIAL  ROTOR  CODE 

The  Full-Potential  Rotor  solves  the  unsteady  full-po¬ 
tential  equation  in  strong  conservation  form.  This  equation 
is  written  below  in  generalized  coordinates. 


±i±\  +  ±t!Z.)  +  LA  +  m  o  (i) 

dr'J*  dV  J  '  dnK  J  '  d(K  J  ’  '  * 


with  density  given  by: 


p={l  +  lyi[-2*r-(U  +  f«)*«- 

(v +  ■».)*,- 1* +  ft)**l}r4r 
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V ,  V  and  W  are  contravariant  velocities  given  by 


U  —  it  +  A  +  A«*,,  +  As*f 

V  =  »jt  +  +  Aj*,  +  A«*f  (3) 

W  =  ft  +  +  Ae4n  +  As#{ 

where 

Ai  =  il  +  £j  +  il 
Aj  =  i>2  +  qj  + 

+  -£  +  4  +  4  (4) 

A«  =  £*»Jx  +  £*»>»  +  i,flM 

Aj  =  £*fx  +  iyfjt  + 

a«  =  v*  u  +  »ivf»  +  n«f« 

All  velocities  are  normalized  by  a^,,  distances  by  the  airfoil 
chord  length,  and  time  by  the  combination  (c/aao).  Density 
is  normalized  by  the  freestream  value. 

Equation  (1)  is  solved  by  using  first-order  backward 
differencing  in  time  and  second-order  central  differencing 
in  space.  The  temporal  density  derivative  is  locally  lin¬ 
earized  about  the  old  time  levels  in  a  manner  that  preserves 
the  conservative  form.  The  resulting  equation  is  approxi¬ 
mately  factored  into  £,  tj ,  and  f  operators.  The  steady-state 
version  of  this  equation  can  be  written  as 


rotor  and  the  attached  finite-difference  grid  move  through 
still  air  and  the  velocity  potential,  ♦,  is  equal  to  zero  in  the 
far  field. 

Boundary  conditions  for  the  calculation  are  shown  in 
Fig.  1.  At  the  outer  radial  boundary  of  the  O-grid,  a  non¬ 
reflection  condition  is  used  to  prevent  the  accumulation  of 
disturbances  there.  This  boundary  condition  is  described 
in  Ref.  9.  On  the  surface  of  the  blade,  a  transpiration  con¬ 
dition  is  used  to  simulate  angle  of  attack  conditions.  Along 
the  inner  boundary  plane,  normal  to  the  rotor,  the  spanwise 
contravariant  velocity  V,  is  set  equal  to  i?t- 

For  lifting  cases,  the  shed  vorticity  is  specified  as  a 
jump  in  potential,  f\  imposed  across  the  coordinate  plane 
i  =  0.  This  coordinate  plane  is  approximately  aligned  with 
the  shear  layer  from  the  trailing  edge.  Density  is  assumed 
to  be  continuous  across  the  cut,  and  an  unsteady  transport 
equation  for  the  potential  jump  is  implemented  in  the  wake. 

The  streamwise  flux  terms  use  upwind  density  bias¬ 
ing  in  regions  of  supercritical  flow  to  ensure  stability  of  the 
algorithm.  The  density  is  biased  in  both  f  and  the  £  di¬ 
rections.  The  density  always  biases  upwind  no  matter  how 
the  grid  lines  are  oriented  to  the  freestream. 
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where  S(,  Sn,  and  if  represent  central-difference  operators 
in  space.  A  steady  state  ADI  relaxation  algorithm  solves 
this  factored  equation.  The  quantity  Rg 0  represents  a  nu¬ 
merical  truncation  error  term  caused  by  incomplete  metric 
cancellation. 

A  sample  grid  for  a  rotor  calculation  is  also  shown  in 
Fig.  1.  An  O-grid  has  been  chosen  for  the  basic  finite- 
difference  grid  because  of  its  efficient  use  of  grid  points. 
Rotor  flows  are  computed  by  assigning  an  appropriate  ro¬ 
tational  coordinate  velocity  to  each  grid  point.  These  coor¬ 
dinate  velocities  are  given  by  £<,  q(,  and  ft.  As  a  result,  the 


If  the  finite-difference  grid  is  extended  radially  be¬ 
yond  the  sonic  cylinder,  then  the  freestream  flow  relative  to 
the  blade  is  supersonic  and  the  full-potential  equation  be¬ 
comes  hyperbolic.  The  FPR  code  uses  central-differencing 
in  space  with  an  upwind-biased  artificial  density  in  regions 
of  supersonic  flow.  This  method  has  been  shown  in  Ref.  11 
to  give  stability  problems  in  regions  where  the  freestream 
flow  is  supersonic. 

In  order  to  remedy  these  instabilities  in  the  far  field, 
the  coefficients  D  and  E  have  been  added  in  the  original 
FPR  algorithm  to  act  as  stabilizers.  The  S  is  a  finite- 
difference  operator  chosen  so  that  it  always  adds  upwind 
data  to  the  system.  The  parameters  D  and  E  control  how 
much  damping  is  added  in  the  £  and  f  directions  respec¬ 
tively.  The  effect  of  these  coefficients  is  to  increase  the 
diagonal  dominance  in  the  solution  matrices.  These  terms 
approach  zero  as  the  solution  converges  to  a  steady-state 
result.  This  modification  to  the  FPR  code  ensures  a  sta¬ 
ble  convergence  at  outer  radial  stations  for  cases  where  the 
freestream  Mach  number  is  well  above  one. 


THE  FAR-F1ELD  GEOMETRICAL 
ACOUSTICS  FORMULATION 

Numerical  data  for  the  potential  in  a  region  near  the 
blade  tip  can  be  used  in  a  geometrical  acoustics  approxima¬ 
tion  for  pressure  in  the  far  field.  This  approximation  will 
not  be  valid  however,  near  the  blade  tip  where  the  flowfield 
is  non-linear  and  rapidly  varying. 

The  geometrical  acoustics  approximation  can  be  based 
on  the  small-disturbance  potential  equation.  One  form  of 
this  equation  is 
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The  time-differentiated  second  order  nonlinearity  Q  is  given 
by 
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The  location  for  the  pressure  evaluation  is  determined 
by  the  values  of  0  and  T ,  where  0  is  a  non-dimensional  ra¬ 
dial  locatation  and  T  is  determined  by  the  aximuthal  angle 
relative  to  the  blade.  The  use  of  Eqn.  (8)  is  restricted  to 
the  plane  of  the  rotor  blade. 

The  surface  integrals  in  Eqn.  (8)  are  located  at  the 
sonic  cylinder,  (r  =  Oq/u).  Choosing  this  location  for  ini¬ 
tial  data  greatly  simplifies  the  result.  The  integrals  at  this 
location  are  evaluated  using  the  computed  results  from  the 
FPR  finite-difference  calculation. 


where  ~i  is  the  ratio  of  specific  heats  and  a0  is  the  sound 
speed  in  undisturbed  air.  Note  that  Eqs.  (6)  and  (7)  assume 
that  the  local  sound  speed  is  constant.  This  is  not  the  case 
however  with  the  finite-difference  solution  close  to  the  blade 
surface. 


The  objective  is  to  use  Eqn.  (6)  to  find  the  disturbance 
pressure  in  the  far  field  in  terms  of  nonlinear  initial  data 
on  some  surface  near  the  blade  tip.  This  is  accomplished 
by  using  a  Kirchhoff  solution  to  Eqn.  (6)  that  converts  it 
to  an  integral  equation  with  Q  as  a  source  term.  Next,  a 
stationary  phase  approximation  is  applied  that  reduces  the 
domain  of  all  integrals  to  that  of  a  data  surface  near  the 
tip.  This  initial  data  surface  is  determined  from  the  finite- 
difference  code  that  solves  the  full-potential  equation  in  the 
region  close  to  the  blade  tip.  The  advantage  of  the  geomet¬ 
rical  acoustics  approximation  is  that  the  finite-difference 
computation  need  not  be  continued  from  the  blade  surface 
all  the  way  into  the  acoustic  field. 


With  these  ideas  in  mind,  Eqn.  (6)  can  be  solved  for 
the  acoustic  pressure  at  a  point  in  the  far  field.  Details 
of  this  solution  are  quite  lengthy  and  will  be  described  in 
Ref.  12.  The  resulting  formula  for  acoustic  pressure  is  given 
by 


INTEGRATION  OF  THE 
FINITE-DIFFERENCE  RESULTS 


The  FPR  solution  produces  pressure  data  at  the  sonic 
cylinder  on  an  O-grid  plane  of  data  perpendicular  to  the 
rotor  axis.  This  solution  is  then  interpolated  onto  a  rect¬ 
angular  mesh  for  ease  of  integration.  Since  this  plane  does 
not  lie  on  a  constant  radial  station,  a  small  approximation 
is  introduced  here  which  greatly  simplifies  the  integration. 
The  integrals  then  use  the  rectangular  mesh  data  for  in¬ 
put.  A  bi-cubic  spline  smooths  the  date  between  points  on 
the  mesh.  A  domain  of  ±3.5  chords  in  both  the  vertical 
and  horisontal  directions  provides  sufficient  input  for  the 
integration. 

All  three  integrals  in  Eqn.  (8)  contain  singularities 
when  T  is  equal  to  o.  In  all  three  cases,  the  equations 
ere  integrated  in  the  o  direction  up  to  a  small  distance  ±< 
from  the  singularity.  The  contributions  at  the  singularity 
are  computed  separately.  The  functions  inside  the  second 
and  third  integrals  can  be  shown  to  be  symmetric  and  odd 
at  the  singularity.  Hence  their  integrated  contributions  are 
*«ro.  The  function  inside  the  first  integral  is  even,  with  a 
non-canceling  contribution  near  the  singularity.  The  value 
of  this  integral  approaches  the  following  limit  at  the  singu¬ 
larity. 
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This  singularity  is  not  insignificant.  It  is  usually  the 
same  order  as  the  rest  of  the  integral  but  of  opposite  sign. 
Near  the  pressure  peak  however,  it  suddenly  switches  sign 
and  reinforces  the  rest  of  the  integral  value.  At  the  pressure 
peak  it  contributes  about  25%  to  the  first  integral. 


Away  from  the  pressure  peak,  the  last  two  integrals  are 
about  the  same  magnitude  as  the  first.  At  the  peak  how¬ 
ever,  the  first  two  integrals  may  double  in  their  respective 
values.  The  third  integral,  however,  increases  its  value  by 
almost  two  orders  of  magnitude,  showing  how  important  it 
is  to  the  overall  solution. 


PREDICTION  OF  IMPULSIVE  ROTOR  NOISE 


Our  rotor-noise  prediction  scheme  is  demonstrated  for 
two  of  the  experimental  test  cases  presented  in  Ref.  13.  In 
this  experiment,  a  1/7  scale  model  of  a  UH-lH  main  rotor 
with  a  NACA  0012  airfoil  section  with  a  root-to-tip  washout 
of  10.9°.  This  model  was  run  at  high-speed  hover  with 
thrust  set  to  a  low  level  in  order  to  minimise  the  wake  in¬ 
fluence.  Acoustic  pressure  data  was  measured  in  the  plane 
of  the  rotor  at  a  radial  location  equal  to  3J2.  A  variety 
of  hover  tip  Mach  numbers  were  tested  from  A#  =  .8  to 
M  =  1.0 

A  top  view  of  the  finite-difference  grid  for  the  M  = 
.9  computation  is  shown  in  Fig.  2.  This  figure  shows  a 
two-dimensional  slice  of  the  grid  in  the  plans  of  the  rotor. 
The  blade  tip  is  located  at  13.7  chords  in  the  spanwise 
direction.  Only  the  outer  three  chords  of  the  blade  are 
computed  in  the  FPR  solution.  A  large  number  of  grid 
planes  are  clustered  in  the  far  field  in  order  to  obtain  good 
resolution  of  the  acoustic  pressures  at  the  sonic  cylinder  (r 
=  15.22  chords).  A  total  of  60  grid  points  were  used  in  the 
chordwise  direction,  38  in  the  spanwise  direction,  and  17  in 
the  direction  normal  to  the  surface  of  the  blade. 


represents  approximately  10  minutes  of  CPU  time  on  a 
Cray  XMP  4/8  computer.  The  addition  of  the  far-field 
damping  terms  causes  this  abnormally  long  run  time.  These 
damping  terms  are  needed  to  control  instabilities  in  regions 
where  the  freestream  is  supersonic.  Typical  iteration  counts 
for  the  FPR  code  on  conventional  grids  are  about  five  times 
lower. 

Computed  Mach  contours  in  the  plane  of  the  rotor  for 
the  M  =  .9  case  are  shown  in  Fig.  3.  The  local  sonic 
Mach  number  contour,  Mi  =  1,  deforms  due  to  the  local 
perturbation  velocity  field  and  touches  the  blade  tip.  The 
strong  local  shock  on  the  tip  portion  of  the  rotor  blade  is 
now  extended  beyond  the  blade  tip  to  the  far  field.  This 
phenomenon  is  called  delocalisation.  Delocalisation  causes 
a  stong  shock  in  far-field  acoustic  signatures  at  tip  Mach 
numbers  beyond  the  critical  delocalisation  Mach  number. 
The  critical  delocalisation  Mach  number  for  this  particular 
blade  geometry  is  0.9. 

Figs.  4a  and  4c  show  computed  pressure  distributions 
from  the  FPR  code  on  the  rotor  plane  at  the  sonic  cylinder 
for  the  M  =  .88  and  M  =  .9  cases  respectively.  Only  the 
in-plane  pressure  distributions  are  shown  in  these  plots  al¬ 
though  the  solution  is  actually  obtained  across  a  plane  at 
the  sonic  cylinder.  These  finite-difference  predictions  are 
then  used  as  initial  data  for  the  Kirchoff  integral  formula¬ 
tion  in  Eq.  (8)  and  the  solution  is  carried  to  the  far  field. 

The  integrations  in  Eq.  (8)  require  initial  data  over  the 
entire  sonic  cylinder.  The  FPR  solutions  only  provide  this 
data  ±3.5  chords  along  the  sonic  cylinder  in  both  the  x  and 
the  s  direction.  Outside  of  this  region,  the  pressures  and 
radial  velocities  are  assumed  to  be  sero.  This  is  clearly  the 
case  in  the  chordwise  direction  as  shown  in  Figs.  4a  and  4c. 
Here  the  pressure  is  seen  to  fall  quickly  to  zero  inside  the 
chordwise  computatinal  domain  for  the  FPR  rode.  Peak 
pressure  along  the  vertical  axis  along  the  sonic  c>  Under  also 
falls  rapidly  to  freestream  values.  This  is  demotistrated  in 
Fig.  5  where  the  relative  pressure  has  dropped  to  zero  well 
within  the  ±3.5  chord  domain  of  the  FPR  code. 


An  untwisted,  nonlifting  NACA  0012  rotor  blade  was 
used  for  the  acoustic  prediction.  This  simplifies  the  finite- 
difference  calculation  because  the  rotor  wake  can  be  ig¬ 
nored.  Boxwell  et  el. 13  have  shown  that  thrust  level  is 
unimportant  for  high-speed  impulsive  noise  cases.  In  gen¬ 
eral  however,  the  FPR  code  can  incorporate  any  of  a  num¬ 
ber  of  rotor  wake  models. 

Approximately  700  iteration  steps  were  necessary  to 
obtain  a  converged  steady-state  solution  for  this  case.  This 


Figs.  4b  and  4d  compare  computed  far-field  pressure 
results  to  experimental  microphone  measurements  at  the 
r  =  3R  radial  station.  Good  agreement  is  seen  between  the 
prediction  and  the  data  for  peak  pressure  magnitudes  and 
pressure  pulse  waveforms.  For  the  M  =  .88  case  (Fig.  4b) 
the  computed  result  underpredicts  the  peak  amplitude  and 
shows  a  slightly  wider  pulse  width  compared  to  the  exper¬ 
imental  data.  The  symmetrical  pulse  shape  agrees  quite 
well  with  the  data  however. 
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The  experimental  data  for  the  M  =  9  cue  (Fig.  4d) 
show  a  strong  shock  in  the  pressure  data  that  is  not  seen 
in  the  prediction.  The  prediction  does  show  a  steeper  slope 
on  the  compression  side  however.  This  discrepency  in  the 
shock  prediction  may  be  attributed  to  two  causes.  First,  the 
Kirchoff  integral  formulation  hu  an  assumption  of  constant 
sound  propagation  speed  as  shown  in  Eq.  (8).  The  far- 
field  nonlinear  phenomenon  may  steepen  the  waveform  as 
described  by  Isom14.  If  this  is  the  case,  then  the  Kirchoff 
integral  formulation  cannot  be  expected  to  model  this. 

A  second  cause  for  disagreement  in  the  shock  structure 
may  be  that  the  FPR  finite-difference  solution  has  smeared 
an  existing  shock  at  the  sonic  cylinder.  Inadequate  grid 
resolution  could  be  the  cause  for  this.  If  this  is  the  case, 
then  the  resulting  far-field  pressure  pulse  will  also  show  a 
smeared  shock.  In  spite  of  these  problems,  the  good  agree¬ 
ment  between  experiment  and  computation  appears  to  val¬ 
idate  our  new  model  for  far-field  rotor  acoustic  prediction 
techniques. 

CONCLUSIONS 

A  new  method  has  been  developed  to  predict  the  im¬ 
pulsive  noise  signature  generated  by  a  high-speed  transonic 
rotor  blade.  The  matching  of  a  finite-difference  method 
and  a  Kirchhoff  integral  formula  combines  the  advantages 
of  each  formulation  in  the  regions  where  each  is  most  ap¬ 
propriate.  The  finite-difference  solution  is  used  close  to  the 
blade,  where  flowfield  nonlinearities  are  large.  The  Kirch- 
hoff  integral  is  used  in  the  far  field  where  its  simplicity 
allows  its  application  over  a  large  domain. 

Predictions  with  this  method  have  shown  excellent 
agreement  with  experimental  data  for  two  hover  cases.  One 
of  these  cases  exhibits  strongly  delocalised  flow  These  re¬ 
sults  represent  the  first  successful  prediction  technique  for 
peak  pressure  amplitudes  with  a  computational  code. 
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Fig.  1  Grid  and  boundary  condition*  for  the  FPR  code. 
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ROTATION  (CHORDS) 

Fig.  2  Top  view  of  the  finite-difference  grid  for  the  M  —  . 9  FPR  calculation. 
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Fig.  3  Predicted  Mach  contours  for  the  M  =  .9  hover  ewe. 
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